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ABSTRACT 

Aims. Our goal is to determine the molecular composition of the circumstellar disk around AB Aurigae (hereafter, AB Aur). AB Aur 
is a prototypical Herbig Ae star and the understanding of its disk chemistry is paramount for understanding the chemical evolution of 
the gas in warm disks. 

Methods. We used the IRAM 30-m telescope to perform a sensitive search for molecular lines in AB Aur as part of the IRAM Large 
program ASAI (A Chemical Survey of Sun-like Star-forming Regions). These data were complemented with interferometric observa¬ 
tions of the HCO^ 1—>0 and C'^O 1—>0 lines using the IRAM Plateau de Bure Interferometer (PdBI). Single-dish and interferometric 
data were used to constrain chemical models. 

Results. Throughout the survey, several lines of CO and its isotopologues, HCO^, H 2 CO, HCN, CN, and CS, were detected. In 
addition, we detected the SO 64—>33 and lines, confirming the previously tentative detection. Compared to other T Tauri and 

Herbig Ae disks, AB Aur presents low HCN 3^2/HCO^ 3—>2 and CN 2—»1/HCN 3—>2 line intensity ratios, similar to other transition 
disks. AB Aur is the only protoplanetary disk detected in SO thus far, and its detection is consistent with interpretation of this disk 
being younger than those associated with T Tauri stars. 

Conclusions. We modeled the line profiles using a chemical model and a radiative transfer 3D code. Our model assumes a flared 
disk in hydrostatic equilibrium. The best agreement with observations was obtained for a disk with a mass of 0.01 Mq, Rin =110 AU, 
Rout =550 AU, a surface density radial index of 1.5, and an inclination of 27°. The intensities and line profiles were reproduced within 
a factor of ~2 for most lines. This agreement is reasonable considering the simplicity of our model that neglects any structure within 
the disk. However, the HCN 3—>2 and CN 2—>1 line intensities were predicted to be more intense by a factor of >10. We discuss 
several scenarios to explain this discrepancy. 

Key words, stars: formation - stars: individual: AB Aur - stars: pre-main sequence - stars: variables: T Tauri, Herbig Ae/Be - 
circumstellar matter - protoplanetary disks 


1. Introduction 

Circumstellar disk s are commonly observed around pre-ma in se¬ 
quence stars (e.g., [Howard et aklbooT: IStrom et aHllOSOl) . The 
formation of disks, together with ejecta phenomena such as out¬ 
flows and jets, dissipate away the excess of angular momentum 
that prevents accretion from the parent cloud. The chemical com¬ 
position of dust and gas contained in these disks provides infor¬ 
mation a bout the initial con ditions in the formation of planetary 
systems (iDutrev et al.ll20l4) . 

The comprehension of chemistry in disks is an important step 
toward understanding of the formation of complex organic, even 
prebiotic molecules on planets. However, the disk chemistry is a 
very unexplored field from the observational point of view with 


very few molecular detections. This scarcity of mol ecules seems 
more accentuated in disks around Herbig Ae stars dOberg et al.l 
I 2 OI ll) . This is mainly due to the low molecular abundances in 
a gas disk that itself has a low mass content. The ultraviolet ra¬ 
diation from the central star photodissociates molecules in the 
surface layers of the disk. Deeper in the midplane, the tempera¬ 
tures drop, and all the detectable molecules freeze out onto dust 
grains. As a result, molecules can only survive in the gas phase 
inside a thin layer. For F and A stars with effective temperatures 
in the range between 6000 to 10000 K, the UV-photons pene¬ 
trate deeper into the disk than the colder M and K stars (Teff 
~ 2500 - 5000K), causing the drop of the molecular detection 
rates. Most species detected are simple molecules, molecular 
radicals, and ions, such as CO, ^^CO, C^^O, CN, CS, C^'^S, C 2 H, 
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H 7 D+, N 7 H+, C-C 3 H 7 , H 2 CO, H 7 O, and HD (e.g.. 

- i - 2 , 

Kastner et al.l 

19971 van Dishoeck et al. 2003: iThi et al. 2004: Oi et alJ 20081: 

Guilloteau et al.ll2006t Pietu et al.ll2007 : Dutrev et 

^12007). 


Unfortunately, most disks remain unresolved even with the 
largest millimeter interferometers. A detailed study of the chem¬ 
ical composition of the gas in a disk requires not only high an¬ 
gular resolution observations in dust continuum and molecular 
lines, but also accurate chemical and physical models. These 
models help to constrain the disk structure in accordance with 
the observations and calculate the molecular abundance pro¬ 
files. In recent years, the chemical and radiativ e transfer mod¬ 
els have improved their performances, (s e e, e.g., Thi_eUa J I2013 


Pinte_e^ai]_l20ld iNomura et alJ 120091: lAgundez et al 
iDutrev et alj|2007h . 


2008 


Herbig Ae/Be stars are intermediate-mass pre-main se¬ 
quence sources that emit much stronger thermal UV radiation 
than do T Tauri stars (TTs). Therefore, their crrcumstellar disks 
are warmer and more ionized. Our target, AB Aurigae (here¬ 
after, AB Aur) is one of the best-studied Herbig Ae stars that 
host a prototypical Her big Ae disk. It has a spectral type AO-Al 
(iHerndndez et al.ll2004l) . It has a M*~ 2.4 Mp , a Teff ~ 9500 K, 
and i t is located at a distance of 145 pc dvan den Ancker et al.l 
Il998l) . The disk around AB Aur shows a complex structure. It 
is larger (Rout~1100 AU, ~7"), than those around TTs, and it 
shows spiral-arm features traced by millimeter contin uum emis¬ 
sion, at about 140 AU from the star (iPietu et al.ll2005h . 

By modeling the '^CO and its isotopologue lines obtained 
from subarcsec imaging with the Plateau de Bure Interferome¬ 
ter (PdBI), it was found that, contrary to typical disks associ¬ 
ated with TTs, the AB Aur disk is warm (>2 5 K all across th e 
disk) and shows no ev idence of CO depletion (iPietu et al.ll2005h . 
ISchrever et alJ (l2008h carried out a chemical study of the disk 
around AB Aur with the PdBI searching for the HCO^ 1—>0, CS 
2—>1, HCN 1—>0, and C 2 H 1—>0 lines, but with only one detec¬ 
tion, the HCO^ 1—>0 line. They propose that the poor molec¬ 
ular content of this disk i s because the UV-p hotons dissociate 
molecules. More recently. Ipuente et alJ (120101) have carried out 
a molecular search using the IRAM 30-m telescope. As a result, 
they detected the SO 34—>23 line. SO had never been detected 
before in a T Tauri or Herbig Ae disk. Its detection provided 
further support to the interpretation of a warmer chemistry in 
this kind of Herbig Ae disk, even if there are some observations 
that suggest the possibility that the SO emission could be com¬ 
ing from an outflo w or from an envelop e rather than from the 
circumstellar disk (iGuilloteau et al.ll2013h . 


In this paper we investigate the molecular composition on 
the AB Aur disk from data obtained as part of the IRAM Large 
program ASAI (A Chemical Survey of Sun-like Star-forming 
Regions) (Pis: R. Bachiller, B. Lefloch). These data were com¬ 
plemented with interferometric observations using the IRAM 
PdBI. Single-dish and interferometric data were used to con¬ 
strain chemical models. A spectral survey is a powerful tool for 
characterizing the chemical composition of a source. It is the 
only way to get a complete census of all molecular species and 
to provide several lines of the same molecule, giving the possi¬ 
bility of multiline analysis and modeling. 

To reach our ai m, we use an updated v ersion of the chemical 
model described in lAgundez et alJ (|2008|) and the 3D radiative 
transfer code MCFOST (see, e.g., ??). 


Table 1. Observed band ranges with IRAM 30-m telescope. 


Band 

Erequency 

[GHz] 

rms 

(mK) 

Av 

(km s“') 

HPBW 

(") 

rib 

E090 

84.5 - 96.3 

5-7 

0.6 

29-25 

0.85 

E150 

133.8- 144.8 

5-8 

0.4 

18-17 

0.79 

E230 

200.5- 272.0 

5-16 

0.3 

12-9 

0.67-0.6 


2. Observations 

2.1. IRAM 30-m radiotelescope 

We carried out a spectral survey toward AB Aur 
(Q'j2ooo=04h55m 45.8 s, (5j2ooo=30° 33’ 04”.2) as part of 
the IRAM Large program ASAI using the IRAM 30-m tele¬ 
scope at Pico Veleta (Granada, Spain). Several observing periods 
were scheduled on July and February 2013 and January and 
Mars 2014, to cover the 1, 2, and 3 mm bands shown in Table[T] 
The Eight Mixer Receivers (EMIR) and the fast Eourier Trans¬ 
form Spectrometers (ETS) with a spectral resolution of 200kHz 
were used for the whole project. The observing procedure was 
wobbler switching with a throw of 200" to ensure flat baselines 
and to avoid possible contamination from the envelope toward 
this young disk. In Table 1 we show the beam efficiency, half 
power beam width (HPBW), spectral resolution, and sensitivity 
achieved in each observed frequency band. 

The data reduction and the line identification were carried 
out with the package CLASS of GILDAS software dMaret et al.l 
I 2 OI ll) . Three databases were used to identify the lines: (1) 
Cologne Database for Molecular Spectroscopy (CDMS; see 
iMuller et al.ll2005l) . (2) the molecul ar spectroscopy da tabase of 
Jet Propulsion Labora tory (J PL; see lPickett et al.lfT998h . and (3) 
MADEX (ICernicharoll2012l) . 

2.2. Plateau de Bure interferometer 

The interferometric observations were carried out in the sec¬ 
ond half of 2011 using the C and D configurations of 
the PdBI with six antennas. This configuration provided a 
beam of a!5.41"x4.31" with a position angle (PA) of 129° 
at the frequency of HCO^ 1—>0 (89.1885 GHz) and of 
a:5.41"x4.31" with a PA of 139° at the frequency of C'^O 1^0 
(112.3593 GHz). During the observations, one 20 MHz band¬ 
width correlator unit was placed at the frequency of the HCO^ 
1^0 line, providing a spectral resolution of ~39 kHz (= 0.13 
km s '). The C'^O line was observed with the WideX correlator, 
which provides a spectral resolution of ~2 MHz (= 5.4 km s“'). 
We used 3C454.3, MWC349, 3C84, 0552H-398, and 0415 -h 379 
as phase and flux calibrators. 

Data reduction and image synthesis were carried out using 
the GILDAS software. The channels free of line emission were 
used to estimate the continuum flux that was subtracted from 
the spectral maps. To improve the S/N, the HCO^ image was 
created with a velocity resolution of 0.4 km s *. The rms of the 
resulting cube was ~6 mJy/beam. In the case of C'^O, we kept 
the original spectral resolution and produced a cube with rms of 
~2 mJy/beam. 


3. Results 

3.1. IRAM 30-m 

In Table ^ we show the list of molecular lines detected at a 
level >5cr. This high S/N is required for a secure detection and 
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Fig. 1. Spectral lines detected toward AB Aur with the IRAM 30-m telescope. The transitions CO 2—>1, HCO^ 1—>0, SO 34 —> 23 , CN 1—>0, 
and CS 3—>2 were taken from iFuente et al.l ll2010h . where cloud velocities were blanked. The HCN 3—>2 is an averaged spectrum from previous 
observations by Fuente et al. (2010) and ASAI data. Magenta line shows the narrow C'*0 1—>0 line profile that traces the emission from the 
foreground cloud, and for visualization purposes the C'*0 1—»0 line (Tpeak~ 0.7 K) is scaled to fit in each panel. 


for having good line profiles for our analysis. For completeness. 



Fig. 2. Comparison of the averaged spectra of the SO and H 2 CO lines 
with the FIC0^3^2 line profile. In the top of these panels we have the 
HC0^3^2 profile showing the double-peak shape that is characteristic 
of a rotating disk. In the middle we can see the averaged H 2 CO profile 
that keep the same double-peak feature in the same range of velocities. 
In the bottom the averaged SO profile is in complete agreement with the 
two other spectra. 


we also include those lines previously detected by Fuente et al. 

( 2010 ). 

The spectra are shown in Fig. [T] The lines of CO and 
'^CO show the narrow line feature (Av ~ 0.5 km s ') typical 
of the lines coming from the foreground cloud (see discussion 
in Fuente et al. 2010). This suggests that their integrated line 
emission is dominated by the cloud emission. However, all the 
other lines show the double peak profile characteristic of rotating 
disks. In Sect. 13.31 we compare the flux of the HCO^ 1^0 line 
listed in Table with the one obtained from our interferometric 
map and demonstrate that the HCO^ 1—»0 emission is coming 
from the disk. For this reason, we adopt HCO^ as a pattern for 
the disk emission. 

The AB Aur disk presents some peculiarities. In a first eye 
verification of the detected lines, we realized that HCN 3^2 
line in AB Aur has an intensity comparable to the less common 
species like H 2 CO. Oberg et al. (2010, 2011) carried out a disk 
imaging survey with SMA toward twelve sources that cover a 
wide range of stellar luminosities among the TT and Herbig Ae 
stars. They show that HCN is strong in disks: HCN 3—>2 and CN 
2^1 are the strongest lines after HCO^ and CO, and are twice 
as intense than the H 2 CO at 1 mm. 

Besides this, SO has never been detected in a protoplanetary 
disk. It is, however, very a bundant in the shocks produced b y 
bipolar outflows, (see, e.g., iBachiller & Perez Gutierrezlll997h . 
In AB Aur disk, we have detected two transitions of SO, 54 ^ 43 , 
and 56 —>45 in the 1mm band, and Fuente et al. (2010) detected 
the SO 34—>23 at 2mm. All the SO lines present double-peak line 
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Table 2. Detections in AB Aur disk. 


Molecule Transition 

Erequency 

[GHz] 

Eup 

[K] 

Area* 

[mKkm 

rms** 

‘ImK] 

'^CO 

2^1 « 

230.538 

16.6 

24388 

35 

'^CO 

2^1 

220.399 

15.9 

3 375 

5 

C'^O 

2^1 

224.714 

16.2 

364 

8 


2^1 

219.560 

15.8 

735 

9 

HCO+ 

1^0 « 

89.188 

4.3 

47 

4 


3^2 

267.557 

25.7 

610 

15 

H 2 CO 

3o3->2o2 

218.222 

21 

75 

7 


3i3—>2i2 

211.211 

32.1 

100 

10 


3i 2—»2ii 

225.697 

33.5 

114 

14 

SO 

54^43 

206.176 

38.6 

29 

5 


56^45 

219.949 

35 

35 

6 


34-»23 0 

138.178 

15.9 

26 

4 

HCN 

3^2 

265.886 

25.5 

144 

14 

CN 

1^0 0 

113.490 

5.5 

26 

6 

CN 

2^1 

226.874 

16.3 

63 

8 

cs 

3^2 0 

146.969 

14.1 

75 

7 


o Transitions previously detected bv iFuente et alj ( l201(]h . * Integrated 
intensity area in the velocity range [3,8] km s“'. ** rms in a channel of 
0.26 km s“' in the 1mm band and 0.65 km s“' in the 3mm band. 


profiles consistent with the interpretation of the emission coming 
from the circumstellar disk of AB Aur (see Fig.ITJ. 

In Fig. |2] we present a comparison between the averaged 
spectra of H 2 CO and SO 1mm lines with the line prohle of 
HCO^ 3—>2. We only averaged the lines observed in the 1mm 
band with similar HPBW, which makes the comparison more re¬ 
liable. The averaged spectra of both H 2 CO and SO match the 
double-peak shape in the same range of velocities as in the case 
of HCO^ line prohle. Waiting for further conhrmation with inter¬ 
ferometric observations, we consider that the circumstellar disk 
is the most likely origin for the SO emission. 

In Table[3j we show the 3 cr upper limits to the column densi¬ 
ties of some interesting species. The column densities were cal¬ 
culated with MADEX (Cernicharo, 2012) assuming a line width 
of Av ~4 km s^^and a source size of ~14 arcsec (2000 AU at 145 
pc). We adopted a molecular hydrogen density of 3.5xl0®cm“^ 
for two representatives temperatures: 30 and 50 K. The hrst one 
is similar to the value derived from the ^^CO observations at 200 
AU by Pietu et al. (2005) and corresponds to the temperature 
in the midplane predicted by our model (see Sect.|5]l. The value 
of 50 K is the mass-averaged kinetic temperature derived from 
our model. We searched for the abundant species HNC, DCN, 
C 2 H, DCO \ N 2 H^, HC 2 N, and C 3 H 2 that have been detected in 
other disks dPutrev et al.ll20l4) . However, we did not detect any 
of these species. They have weak lines whose emission could 
be below our sensitiv ity limit. The tentative detection of H 2 S 
bv iFuente et al.l (1201 Ol) was not conhrmed. We also searched for 
HCS^, HCO, SiO, and CH 3 OH lines, which are abundant species 
in young protostars. They have not been detected coming from 
any disk or around AB Aur. 


ported by iRodrfguez et all (12014h (direction angle~70°). More¬ 
over, the rotation axis of HCO^ is co nsistent with that of the CO 
gas disk traced bv iPietu et al.l (l2005l) (~50°). 

The bottom panels of Fig. [3 show the real part of the visibil¬ 
ities of the AB Aur molecular emission as a function of the uv 


radius. Each black point is an azimuthal average of the visibili¬ 
ties inside circular rings in steps of 5 m from 0 to 200 m. Eor the 
uv analysis, we tried to fit different models, but the data quality 
does not allow us to distinguish between circular and ring disk 
models. Red points indicate the best fit to the uv plot, using a 
circular disk model. In the case of the uv plot of the C'^O emis¬ 
sion, we note that it is well htted by a circular disk model cen¬ 
tered on aj2ooo=04h55m 45.81s, dj2ooo=30° 33’ 03”.83, with 


a flux density of 18+3 mJy and 6±0.8"of diameter. However, 
for the case of the HCO^, the fit is not good at long baselines. 
Although there is much noise in the HCO^ data, the uv plot sug¬ 
gests that there might be two components, a compact unresolved 
component, which is responsible for the flattening of the uv vis¬ 
ibility amplitude at large (>100 m) baselines, and an extended 
component (resolved) with a decreasing trend of the uv visibil¬ 
ity amplitudes with radius at small uv radii. The flux density for 
this HCO^ model was 62+3 mJy, and it had the same diameter of 
5.8+0.3". The fluctuations observed at large mv radii indicate that 
the morphology of the source is complex, consistently with the 
complexity of the structures rev ealed by previous observations of 


l2008tlRodrfsuez et al.ll2014ll. 

Table 3. Non-detections in AB Aur disk.' 

Species 

Erequency 

(GHz) 

rms 

(mK) 

<5v 

(km s“ 

3(T 

')(mK 
km s“' 

N(50K) 

cm -2 

) 

N(30K) 

cm-" 

HNC 

90.664 

4.20 

0.65 

20.32 

2 . 0 x 10 " 

1 . 6 x 10 " 

DCN* 

144.827 

9.50 

0.40 

36.05 

7.0x10'“ 

6 . 0 x 10 '“ 


217.238 

8.00 

0.27 

24.94 

2.3x10'“ 

2.5x10'“ 

CCH 

87.316 

6.20 

0.67 

30.45 

1.9x10'3 

1.4x10'" 

DCO+ 

144.077 

9.20 

0.40 

34.91 

6 . 0 x 10 '“ 

4.6x10'“ 

NjH* 

93.171 

6.50 

0.63 

30.96 

8 . 0 x 10 '" 

6 . 0 x 10 '" 

HC 3 N 

90.979 

4.20 

0.64 

20.16 

3.9x10" 

2.9x10" 


136.464 

7.20 

0.43 

28.33 

1.7x10" 

2.4x10" 

C-C 3 H 2 

217.822 

8.00 

0.27 

24.94 

1 . 0 x 10 " 

1 . 2 x 10 " 


251.311 

7.20 

0.23 

20.72 

8 . 8 x 10 '“ 

1 . 2 x 10 " 

p-HjS 

216.710 

8.60 

0.27 

26.81 

5.4x10'" 

2.7x10'" 

0 -H 2 S 

204.140 

15.00 

0.28 

47.62 

1 . 6 x 10 '“ 

5.2x10'“ 

HCS+ 

85.348 

5.60 

0.69 

27.91 

1 . 8 x 10 '" 

1 . 2 x 10 '" 


256.027 

10.00 

0.23 

28.77 

9.7x10'“ 

1 . 1 x 10 " 

HCOo 

255.341 

10.00 

0.23 

28.77 

1 . 0 x 10 '" 

1.4x10'" 


265.348 

16.00 

0.22 

45.03 

3.2x10'" 

4.5x10'" 

SiO 

86.846 

6.00 

0.67 

29.47 

4.1x10" 

3.2x10'“ 


217.104 

8.50 

0.27 

26.50 

3.5x10'“ 

4.9x10'“ 

CH 3 OH 

252.485 

11.00 

0.23 

31.65 

2.4x10'" 

1.4x10'^' 


Notes. ' The detection limit was calculated assuming a line width of 
Av ~4 km s“'. The upper limits on the column densities was calculated 
with MADEX for a source size of ~14 arcsec (2000 AU at 145 pc), a 
H 2 density of 3.5xl0*’cm“^ for two temperatures: 50 and 30K. * We 
used the collisional rates of HCN, o Assuming LTE. 


3.2. Plateau de Bure 

The upper panels of Eig. [3 show the contours of the integrated- 
velocity intensity maps of HCO^ 1—»0 (on the left) and C'^O 
1—»0 (on the right). Both maps are integrated in the velocity 
range between [4,8] km s“'. The emission coming from the disk 
traced by HCO^ 1 —>0 line is elongated and almost perpendicu¬ 
lar to the emission traced by the collimated outflow recently re- 


3.3. Comparison between IRAM 30-m and Plateau de Bure 

In Eig. |4]we show the HCO^ 1—»0 and C'^O 1—»0 lines. The 
availability of PdB maps allows a comparison between interfer¬ 
ometric and single-dish data. To make this comparison, the in¬ 
terferometric maps have been convolved to the single-dish beam 
size of the corresponding transition. Eor the C'^O 1—>0 line, the 
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Fig. 3. {Upperpanels) HCO^ 1—»0 (left) and C'^O 1—>0 (right) integrated velocity maps obtained with the PdBI from 4 to 8 km s“'. Blue contours 
levels are percentages of the maximum (0.15 Jy beam“‘ for C'^O and 0.22 Jy bearn^* for HCO^) from 10 to 90 in steps of 10. The positi on of 
the central source is marked by a yellow star. The magenta arrow indicates the direction of the jet (~70°) reported bv lRodrfguez et al]ll2014h . The 
red contour represents the half maximum power contour. The ellipse in the left bottom comer indicates the synthesized beam sizes: 5.41" x 4.31" 
with a position angle (PA) of 39° for HCO^ 1—»0 and 4.58" x 3.98" with a PA of 49° for the C*’0. {Lowerpanels) Corresponding uv plots to the 
upper maps. The uv plot shows the real part of visibilities as a function of the uv distance, binned in circular rings in steps of 5m. 


temperature measured with the PdB is about four times less in¬ 
tense than obtained with the IRAM-30 m telescope in the 2^ 1 
line, which is consistent with optically thin emission. For the 
HCO^ 1—>0 line, all the flux measured with the IRAM 30-m is 
recovered by the PdBI, proving that the detected emission comes 
from the circumstellar disk. 

Looking at the HCO^ 1 —>0 spectra of Fig.|4] we note that the 
line profile is doubled peaked, which is characteristic of emis¬ 
sion from a rotating disk. However, the line profile seen by the 
interferometer seems to peak at a slightly blue-shifted velocity 
with respect to the IRAM 30-m spectrum. We believe that this is 
due to the poor S/N in the PdBI spectrum, which means that the 
peak velocity and the other line parameters are affected by large 
uncertainties. 

4. Comparison with other sources 

For a complete picture of the chemistry in the AB Aur transi- 
tio n disk, we compared it w ith the sample of 12 disks observed 
bv lOberg et alJ (l20Toll201 ih . Such disks were entirely observed 
with the SMA telescope, ensuring the uniformity of observa¬ 
tional conditions in the sample. Moreover, the li nes detected 
toward AB A ur are among the lines detected by lOberg et^ 
(1201(1 120111) at 1mm in the other sources, which facilitated 
the comparison. The sample covers a wide range of stellar lu¬ 
minosities among the TT and Herbig Ae stars. It consists of 
six disks in Taurus and six in the southern sky. Among these 
disks, we have several (pre-) transition disks with internal gaps: 
DM Tau, LkCa 15, GM Aur, V4046 Sgr, HD 142527, and SAG 


206462. Table|4]lists the disk types, pT or T for pre- or transition 
disk respectively, and F for full disks. The CN 2—>1/HCN 3—>2, 



V (km/s) V (km/s) 

Fig. 4. Comparison between single-dish IRAM 30-m spectra (upper 
panels) and Plateau de Bure interferometric spectra (lower panels). To 
compare properly both sets of spectra, PdB spectra have been convolved 
to the single-dish beam size of the corresponding transition. For the 
HCO+ 1—»0, the HPBW of the convolving gaussian is 27.0" x 27.2" 
and for the C'^O 1—>0, 9.4" x 9.2" (to obtain the same angular resolu¬ 
tion of the 30-m C*^0 2—>1 line). The conversion factor between Jy and 
K is 7K = 1 mjy/beam for HCO^ 1—»0 and 5K = 1 mjy/beam for C'^O 
1 ^ 0 . 
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HCN 3^2/HCO+3^2, and the H 2 CO 3o3^2o2/HCO+ 3^2 in¬ 
tegrated intensity line ratios are listed in the other columns. Car¬ 
rying out a comparison between several sources is not direct, 
since there are many uncertainties because the position angles 
and inclinations affect which parts of the disk are observed, 
hence our picture of the chemistry. Disk masses are poorly con¬ 
strained, since they are derived from dust emission, assuming 
the canonical interstellar dust/gas ratio value of 1/100. The ac¬ 
tual dust/gas ratio may be variable among the sources and also 
be affected by coagulation and photoevaporation. The difference 
in disk masses may affect the relative abundances of different 
species. 

Among the molecular species detected in AB Aur, we have 
HCN and CN. They are particularly interesting because they 
trace different physical conditions in the photodissociation re¬ 
gions (PDRs) and because their abundances are closely related 
HCN may be ea sily photodissociated by UV radiation from the 
star t o yield CN (Fuente et al.ll993t iBachiller & Perez Gutierre^ 
ll997t[B5ger & Sternberg|l2005h. 

Both molecules, CN and HCN, share the same chemistry 
and are expected to come from the same region, although the 
CN/HCN ratio can yary (see Fig. lA.lb . Furthermore, the CN 2-1 
and HCN 3-2 lines haye similar critical densities and therefore 
trace regions with similar physical conditions. Assuming opti¬ 
cally thin emission, their integrated line intensities are propor¬ 
tional to the number of molecules of each species. With all these 
assumptions, the CN 2-1/HCN 3-2 line ratio is a good tracer of 
the CN/HCN abundance ratio. 

In Fig. |5] we show the integrated intensity line ratios 
and upper and lower limits for CN 2—>1/HCN 3—>2, HCN 
3^2/HCO+3^2, and H 2 CO 3o3^2o2/HCO+ 3^2 ratios. The 
CN/HCN ratios yary by factor 4 between the DM Tau and AB 
Aur. AA Tau is the only signihcant outlier in terms of CN/HCN 
flux ratio, which has a yalue of 2.9. The ayeraged yalue of the 
CN/HCN ratio in the sample is ~1.1. This yalue is similar to the 
one measured in the Orion Bar PDR ~0.85 (J. R. Goicoechea, 
priyate communication) (see Fig. 0. CN and HCN are more 
abundant in the illuminated disk surface (see Fig. lA.ll l. Chemi¬ 
cal models show t hat the CN/HCN ratio is a good tracer of th e 
incident UV held (iFuente et alJll993llBoger & Sternberg|2005h . 
taking yalues of about 1 in dark clouds and increasing to > 3 in 
PDRs. The similarity between the CN 2-1/HCN 3-2 line ratio in 
AB Aur and Orion is consistent with the interpretation of CN 
and HCN emission coming from the PDR formed in the inner 
part of the disk and the illuminated disk surface. 

We found a weak trend between CN 2^1/HCN 3^2 ratio 
and the morphology of the disk. Transition disks with high lu¬ 
minosities haye the lowest yalues for CN/HCN ratios. AB Aur 
is located with the other transition disks, SAO 206462 and HD 
142527, in the bottom right hand part of the plot (see Fig0. 

The ayerage CN/HCN abundance ratio is expected to depend 
on the total area of the illuminated disk surface and the inten¬ 
sity of the UV field. Assuming the same morphology for all the 
disks, one would expect that the CN/HCN ratio increases with 
the stellar luminosity. Howeyer, we hnd the contrary trend. We 
interpret that this trend might be due to the different disk mor¬ 
phologies. In the transition disks with large inner gaps, the emis¬ 
sion is dominated by the outer ring, i.e., the region farthest from 
the central star, where the FUV flux is weaker. This effect would 
be increased if the disks were flatter, because the illuminated 
surface would be drastically reduced. 

In Fig. |5] we also plot the HCN 3^2/HCO+3^2 and 
the H 2 CO 3o3—>2o 2/HCO^ 3—>2 ratio. The ayeraged yalue 
of HCN/HCO^ ratio in the sample is 0.43 and for H 2 CO 


3o3—>2o 2/HCO^ 3—>2 is 0.1. We did not find any trend for 
each of these ratios with the stellar luminosity or the disk mor¬ 
phology. Howeyer, we found a relation between them both. 
Clearly, the HCN 3—>2/HC0^3—>2 ratio is higher than H 2 CO 
3o3—>2o 2/HCO^ 3^2 ratio in all the cases. This is also true 
for the Orion Bar: their corresponding yalues are 0.74 for the 
HCN/HCO+ and 0.11 for H 2 CO/HCO+ (J.R. Goicoechea, pri¬ 
yate communication). 


5. Modelization 


5.1. Physical model 

For deeper insight into the chemistry of this Ae disk, we car¬ 
ried out seyeral modeling tasks. Our model assumes a flared 
disk with the density prohle in the yertical axis determined by 
the assumption of hydrostatic equilibrium. The density and tem¬ 
perature distributions are shown in Fig. IA.2I The temperature 
and yisual extinction are com puted using the RADMC code 
(iDullemond & Dominildl2()()4ll . which solyes the continuum ra- 
diatiye transfer throughout the dis k. We use a dust model with 
the M RN grain size distribution dMathis, Rumnl. & Nordsiecfl 
Il977h made of 100% silicates. Gas kinetic temperature is higher 
than dust temperature in the yery external layers on the disk sur¬ 
face, where the photoelectric effect is the main heating mecha¬ 
nism. Howeyer, the millimeter lines studied in this paper come 
from a denser and cooler region where gas and dust are expected 
to be thermally coupled. In our model we assumed that the gas 
and dust temperature are equal. 

_ Flared disks are expected around young stars (iMeeus et al.l 

l2001h . The input parameters are the disk mass (Mdisk), the in¬ 
ner and outer radii (Rin and Rout) and the radial index, a, of the 
surface density profile (£(R) oc R^“). 

As an initial trial, we considered Models A, B, C, and D 
with the parameters shown in Table 0 We adopted the disk 
mass of 0.02 Mq deriyed by iPietu et al.N2005h from continuum 
1mm PdB obseryations and allowed it to change by a factor of 
2, so within its uncertainty. The yalues of Rin and Rout are yar- 
ied between the yalues found from '^CO a nd those obtained 
from the rare CO ip topologues and HCO^ (iPietu et akl 120051 : 
ISchreyer et al.ll20dM this paper). The yalue of a was taken 2.15 , 
in a greement with preyiou s measurements by IPietu et al.l (l2005h 


andy 

The parameters adopted in Models A,B,C, and D are based 
on fitting interferometric obseryations with a simple disk model. 
More recent high spatial resolution obseryations of CO and '^CO 
showed that the AB Aur disk has a complex structure, with a ro¬ 
tating ring that extends to 4" (560 AU) and a spiral arm e xtend- 
ing up to 7" (^1000 AU) (iTang et al.lUOH [Pantin et al.l[200l . 
Moreoyer, according to obseryations by lHashimoto et aUdMll, 
there are two dusty rings: one inmate compact ring (R « 40 AU) 
with an inclination angle of approximately 40° relatiye to the 


Table 5. Input parameters used for modeling the AB Aur disk. We var¬ 
ied the disk mass, the inner radius, the outer radius of the disk, and the 
radial index, a, of the surface density profile, S(R) oc The disk 
mass is given in units of solar masses (Mg ) and radii in astronomical 
units (AU). 



A 

B 

C 

D 

E 

Mdisk 

0.02 

0.01 

0.01 

0.01 

0.01 

Rin 

70 

70 

110 

70 

no 

I^OUt 

1100 

550 

550 

1100 

550 

a 

2.15 

2.15 

2.15 

2.15 

1.5 
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Table 4. Disk list of the sample of lOberg et al.ll201(il201lL and AB Aur with their corresponding integrated line intensity ratios of the transitions; 
CN 2^1/HCN 3^2, HCN 3^2/HCO+ 3^2, and HjCO 3o3^2o2/HCO+ 3^2. 



Disk 

L. 

CN 2^1 

HCN 3-^2 

H7C03o3^2o2 

Sources 

type 

Lq 

HCN 3^2 error 

HCO+ 3^2 error 

HCO+ 3->2 error Refs. 


DM Tau 

PT 

0.3 

1.62 

0.09 

0.55 

0.03 

0.07 

0.01 

a 

AA Tau 

F 

0.59 

2.98 

0.81 

0.25 

0.08 

0.07 

0.04 

b 

LkCa 15 

pT 

1.2 

1.23 

0.07 

1.06 

0.08 

0.13 

0.02 

a 

GM Aur 

pT 

0.93 

0.76 

0.13 

0.25 

0.03 

0.11 

0.02 

a 

CQ Tau 

F 

8 

>0.44 


<0.67 


<0.28 


c 

MWC 480 

F 

11.5 

1.41 

0.17 

0.50 

0.07 

<0.06 


d 

IMLup 

F 

1.3 

1.01 

0.11 

0.36 

0.04 

0.05 

0.01 

e 

AS 205 

F 

4 

- 


<0.14 


<0.06 


f,g 

AS 209 

F 

1.5 

1.65 

0.18 

0.50 

0.05 

<0.05 


g,h 

V4046 Sgr 

pT 

0.5 

1.22 

0.05 

0.87 

0.03 

0.09 

0.02 

i 

SAO 206462 

T 

8 

1.00 

0.50 

0.49 

0.22 

<0.25 


j 

HD 142527 

T 

69 

- 


<0.10 


0.03 

0.02 

j 

AB Aur 

T 

38 

0.44 

0.10 

0.24 

0.02 

0.12 

0.01 

k 


Notes. (a) lEspaillat et al.| j2010|) ; (b) [White & GhezI ( 1200 ll) ; (c) ]Mannings & Sargentid 19971) ; (d) l^imon et ahll 20001) 7(6) iHughes et al.l lll994h : (f) 
lPratoetalJll200.3h : fgi I Andrews etaP ll2009li : thl lHerbig & Beiil(ll988h : Cii lOuast et all ll2000h : til lGarcia-Lopez et alJll2006h : tkl lDeWarf etaP 
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Fig. 5. Integrated line intensity ratios and upper and lower limits for CN 2—>1/HCN 3—>2 (left) and HCN 3—>2/HCO^ 3—>2 and H 2 CO 
3o3—»2o2/HCO^ 3—>2 (right) transitions, as a function of the stellar luminosity. Horizontal dashed lines indicated the corresponding values for 
the Orion Bar (J. R. Goicoechea, private communication). The letters correspond to each of the sources. (A) DM Tau; (B) V4046 Sgr; (C) AA Tau; 
(D) GM Aur; (E) LkCa 15; (F) IM Lup; (G) AS 209; (H) AS 205; (I) SAO 206462; (J) CQ Tau; (K) MWC 480; (L) AB Aur; (M) HD 142527. 


plane of the sky and th e outer ring with R » 100 AU and an 
inclination angle of 27°. iHonda et al.l (l2010h derived a value of 
a close to 1 for the outer ring. Our model is too simple to ac¬ 
count for this complex structure, but we added Model E (see 
Ta ble Ell to account fo r the lower radial spectral index measured 
bv lHo^a et al.l (l2010l) . 

In all the models, the stellar parameter s were kept fixed to T ,. 
= 10 000 K, R. = 2.5 Ro, M. = 2.4 Mq (ISchrever et al.ll2008h . 
and the gas-to-dust mass ratio to 100. 

5.2. Chemical model 

We used a simple time-dependent chemical model in which the 
representative chemical abundances evolve under f i xed p hys¬ 
ical conditions (u pdated version of lAgilndez et alJ l2008l and 
iFuente et al.ll2010h . The chemical model includes the elements 
H, C, N, O, and S. As i nitial abundances w e adopt the so- 
called “low metal” values (IWiebe et al.l 120031) . with all carbon 
in the form of CO, the remaining oxygen in the form of wa¬ 


ter ice, nitrogen as N 2 , sulfur as CS, and a ionization fraction 
of lxl0~^(typical of dark c louds, see, e.g.. ICaselli et al.l[i998t 
lAgilndez & Wakelamll2013b . Th e gas phase chemical n etwork is 
an extension of the one used by lAgilndez et al.l (l2008h . We also 
consider adsorption onto dust grains and desorption processes, 
such as thermal evaporation, ph otodesorption, and desorpti on in¬ 
duced by cosmic rays (see, e.g.. [Hasegawa & Herbsll[l99^ . The 
photodissociation rates of H 2 and CO are computed including 
self-shielding effect s. For the stellar FUV field we adopt the ap¬ 
proach described in lAgilndez et al.l (l2008l) . while in the case of 
the interst ellar FUV field, w e use the shielding coefficients cal¬ 
culated by iLee et al.l (Il996l) . The adopted H 2 formation rate on 
the grain surface is 3x10^'^ cm^ s“'. 

Using the chemical composition computed by our model, 
we performed the 3D radiative transfer code MCFOST (??) to 
compute the synthetic molecular emission maps. MCFOST uses 
the chemical abundances and kinetic temperature to compute the 
level populations using a Monte Carlo met hod. We followed a n 
algorithm similar to the one published by iHogerheiidel (l2010l) . 
separating the contribution of the local and external radiation 
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fields in each cell. The synthetic emission maps and spectra are 
then produced by a ray-tracing algorithm that formally integrate 
the source function computed by the Monte Carlo method. We 
iterated this procedure several times over a range of position an¬ 
gles from 0° to 40°, relative to the plane of the sky. Finally, we 
convolved the resulting map with the beam of Pico de Veleta to 
compare with observations. The results are shown in Fig. | 6 ] and 
in Figs. lA.3l and lA.4l 

5.3. Comparison between model and observations 

Our models from A to E are reasonable approximations of the 
AB Aur disk (see Figs.l6l lA.3l and lA~4l i. Model A corresponds 
to a disk with mass of 0.02 Mq and a large outer radius of 1100 
AU. This model overestimates the observed intensities and the 
line widths for all transitions. In Model B we used the same pa¬ 
rameters as in Model A, but we changed the outer radius to 550 
AU and the disk mass to 0.01 Mq . By decreasing and redis¬ 
tributing the mass in a smaller radius, we observed a decrease 
in the line intensities, but there is no change in the width of the 
profiles. In Model C we increased the size of internal radius, the 
synthesized spectra fit the width of the lines best, and the line 
profiles are more similar to those of the observed spectra. Still, 
under the assumption of a mass of 0.01 Mq and when testing 
the same other parameters as in Model A, we noted the effect of 
decreasing the mass in the disk in Model D. The lines are still 
wider than observed, although the intensities are adjusted rela¬ 
tively well within a factor of 2 for most of the transitions. 

Finally, Model E is like the previous Model C, but it changes 
the exponent of the radial surface density law from 2.15 to 1.5. 
The best results were obtained for this model, which corresponds 
to the values Mdisk=0.01MQ , Rin =110 AU, Rom =550 AU, and 
a surface density radial index q;= 1.5 with an inclination angle 
of i=27°. Almost all modeled lines agree with the observations 
within a factor of 2 , and the modeled profiles successfully repro¬ 
duce the double peak profiles from the disk emission. The excep¬ 
tions are the HCN 3^2 that is predicted with an intensity around 
30 times higher than the observed ones and the CN 2—»1 with a 
predicted intensity that is 17 times higher. This large difference 
was highlighted in the eye inspection during the data analysis. 
We noted that HCN 3—>2 and CN 2—>1 were especially weak. 
This peculiarity in the emission is also observed as a common 
trend in transition disks. 



Fig. 7. Vertical column densities of several molecules as a function of 
radius at 2.5 Myr as calculated for the chemical model E. 


Table 6. Initial abundances in Models MO, Ml, and M2'. 


Species 

MO 

abun/Htot 

Ml 

abun/Htot 

M2 

abun/Htot 

gas phase 

H 

O.OOOE-hOO 

2.1414E-04 

1.6066E-04 

H 2 

0.500E-H00 

4.9989E-01 

4.9992E-01 

He 

0.085E-H00 

8.5000E-02 

8.5000E-02 

C 

O.OOOE-hOO 

1.2132E-05 

4.5720E-08 

O 

O.OOOE-hOO 

1.1569E-04 

8.6790E-06 

N 

O.OOOE-hOO 

2.1795E-05 

5.1383E-07 

S 

O.OOOE-hOO 

8.6586E-08 

1.5544E-08 

OH 

O.OOOE-hOO 

5.3736E-08 

4.3489E-07 

H 2 O 

O.OOOE-hOO 

2.7982E-07 

1.3307E-07 

O 2 

O.OOOE-hOO 

1.1912E-07 

6.3504E-06 

CO 

7.860E-05 

5.9802E-05 

4.3019E-05 

CO 2 

O.OOOE-hOO 

1.0848E-08 

6.9679E-08 

H 2 CO 

O.OOOE-hOO 

1.0260E-08 

8.1825E-10 

NH 3 

O.OOOE-hOO 

2.6452E-09 

4.6800E-08 

N 2 

1.235E-05 

1.3839E-06 

8.6479E-06 

CN 

O.OOOE-hOO 

1.4022E-08 

4.8573E-09 

HCN 

O.OOOE-hOO 

1.2226E-08 

6.2810E-09 

HNC 

O.OOOE-hOO 

6.988 lE-09 

5.4021E-09 

CS 

9.140E-08 

2.3049E-09 

2.3040E-10 

SO 

O.OOOE-hOO 

3.3169E-10 

1.0105E-08 

ice 

CH 4 

O.OOOE-hOO 

2.3428E-08 

6.1490E-08 

H 20 

1.014E-04 

6.7127E-08 

1.0762E-06 

02 

O.OOOE-hOO 

1.0774E-08 

3.8826E-05 

CO 

O.OOOE-hOO 

3.7799E-06 

3.5103E-05 

C 02 

O.OOOE-hOO 

2.2142E-09 

2.1759E-07 

H 2 CO 

O.OOOE-hOO 

2.6391E-09 

5.4471E-09 

NH 3 

O.OOOE-hOO 

4.8238E-10 

1.4827E-07 

N 2 

O.OOOE-hOO 

3.7448E-08 

3.0605E-06 

'Notation: 2.1414E-04: 

= 2.1414x10^“^ 

We call that our chemical model does not consider surface 


reactions that could affect the HCN and CN chemistry. However, 
our model gives vertical HCN column densities of ~ 10'"^ cm“^, 
in agreement with other more complex mo dels that include sur¬ 
face chemistry (see, e.g. JWalsh et al.ll2010h . Therefore, we think 
that the discrepancies between observations and our predictions 
are more likely due to the assumed morphology of the AB Aur 
disk. 

In Eig.|7]we show the HCN and CN vertical column densities 
as a function of the disk radius. Both species have peaks close 
to the inner radius of 110 AU. Subsequently, they decrease to 
a minimum value around 200 AU and increase again for larger 
radii because of the PDR formed on the surface of the flared 
disk. In the case of a flat disk, the HCN and CN column densities 
would be lower for large radii and in better agreement with our 
observations. This could also be the cause of the low HCN and 
CN intensities in other transition disks. 

In our model the SO synthesized lines agree with the obser¬ 
vations by a factor of 2. SO emission has been dete cted asso¬ 
ciated to disks of young low-mass Class I protostars (lYen et al.l 
|2014|) . The youth of the AB Aur disk could be key for a higher 
SO abundance. Moreover, the AB Aur disk is warmer than those 
associated with low-mass stars and, unlike to these disks, which 
present a standard abundance of CO without s igns of heavy de¬ 
pletion (iPietu et al.ll2005l : ISchrever et al.ll2008h . 

In general, there is a discrepancy of almost a constant factor 
of 2 between our model predictions and observed line intensi- 
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Fig. 6. Comparison between modeled spectra and those detected by 30m telescope toward AB Aur disk. In colors, we have the lines profiles that 
correspond to the synthetic line fluxes performed by MCFOST code. The blue and magenta lines correspond to the same m odel with inclinatio n 
angles for the disk of 27° and 40°, respectively. The vertical dash line indicates the ambient cloud velocity, Visr = 5.9 km s“' dDuvert et alJl986h . 
All the observational spectra show the double-peak shape that is characteristic of a rotating disk. All modeled lines agree with the observations 
within a factor of 2. HCN 3—>2 and CN 2—>1 present the largest differences between the predicted intensities and the observations, 30 and 17 times 
larger than observations, respectively. 


ties. This factor may be due to a different gas/dust ratio from the 
standard value, and it could even vary within the disk. This could 
produce a decrease in the intensity of the lines coming from the 
gas-deficient region. It is definitely not the unique explanation, 
but taking the large uncertainties of our simple model into ac¬ 
count, we think it is not worthwhile discussing this disagreement 
further. We call that the gas/dust ratio is one of the least known 
parameters in all the disk models. 

The shortcoming of time-dependent models is that the pre¬ 
dictions depend on the initial conditions. In the following sec¬ 
tion, we discuss the impact of such initial conditions on our re¬ 
sults. 


5.4. Dependence on initial conditions 

To analyze the impact of the initial condition on our model, we 
used a trial of initial abundances: MO (the actual model used 
up to now in this work). Ml and M2. The values for each set 
of abundances are shown in Table 0 Model MO corresponds to 
a standard solar elemental composition corrected with a deple¬ 


tion around 100 for Na, Si, Mg, and Fe (see lWiebe et al.ll2003l: 
ISemenov et al.ll2004ll . In this model, all the N is in the form of 
N 2 , the C in CO, the O in water ice, and all the S in CS. (We 
chose this option as the limit case to verify the model’s ability 
to form SO.) Ml and M2 initial abundances values correspond 
to the values of a dark cloud (nH=2x lO^cm^^, T(K)=10K) at a 
time of 10^ and 10® years, respectively. 

Figure H] shows the comparison between the observed spec¬ 
tra and the profiles synthesized with the radiative transfer code 
MCFOST for models MO, Ml, and M2 and for an angle of 27°. 
Models Ml and M2 show practically the same line profiles for 
all transitions. As an outstanding result, the CS and SO predicted 
line profiles do not vary with the initial conditions, supporting 
the robustness of our results. 

In the case of CN and HCN, the predicted intensities are 
lower than in MO, but remain more intense than the observa¬ 
tions. Thus, varying the initial conditions does not improve our 
results. This discrepancy between models and observations may 
be due to a different physical structure (see Sect. 15.31 1 or to the 
uncertainties in the nitrogen chemistry. For instance, the values 
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Fig. 8. Detected spectra toward AB Aur overlaid with the synthesized transition lines obtained with the radiative transfer code MCFOST for 
models MO, Ml, and M2 for an angle of 27° for the physical parameters of model E. 


of elemental abundances, the ortho-to-para H 2 ratio, the binding 
ener gies, and time-dependent effects could affect ou r results (see, 
e.g./ lLe Gal et al.l 120141: l^hwarz & Berginll2014h . However, a 
detailed study of these effects is beyond the scope of this work. 

Models Ml and M2 fail to predict the abundances of 
formaldehyde by a factor of 10, because the formaldehyde emis¬ 
sion mainly comes from the midplane where the chemistry pro¬ 
cesses are slow and the molecular abundances are far from the 
equilibrium. In Models Ml and M2, there is, hrst, a period of 
very early formation of molecules, followed by a time of destruc¬ 
tion of the chemical complexity (between a few 10^ and 10^ yr), 
and then the formaldehyde abundance increases again to reach 
the equilibrium value at 10^ yr. The steady-state abundance is 
more like the result of our MO model. 

Among the studied molecules, only the abundance of H 2 CO 
strongly depends on the assumed initial conditions. 

6. Summary and conclusions 

We present in this paper the molecular survey toward the proto¬ 
typical disk around AB Aur. Our results and conclusions can be 
summarized as follow. 

1. We detected the 1—>0 and 2—>1 rotational transitions of 
'^CO and its isotopologues: '^CO, C'^O, and C^^O. We also 


detected the lines HCO^ 1—>0 and 2^1, H 2 CO Sqs—> 2 o 2 , 
3i3^2i 2, and 3i2^2ii, SO 54^43, 56^45, CS 3^2, HCN 
1—>0, 3^2, and CN 1—>0 and 2—>1. Based on the line pro- 
hles and on the agreement between the velocity ranges, we 
interpreted that the emission from all these lines comes from 
the disk. We therefore conhrmed the previous detection of 
SO in AB Aur by Fuente et al. (2010). 

2. Some species were not detected in our survey, even if they 
were observed in other disks, such as HNC, DON, C2H, 
DCO^, N2H^, HC3N, and C3H2. In general, these species 
have weaker lines whose emission could be below our sensi¬ 
tivity limit. 

3. Carrying ou t a comparison with other 12 disks (lOberg et al.l 
I 2 OIOII 2 OI ill that cover a broad luminosity range and disk 
evolutionary stages, we found that AB Aur presents the typ¬ 
ical chemical signature of a transition disk. Transition disks 
have lower CN 2—>1/HCN 3—>2 ratios than full disks. More¬ 
over, the HCN 3^2/HCO+ 3^2 and H 2 CO 3^2/HCO+ 
3^2 ratios are similar. 

We modeled the line prohles using the chemical code by 
lAgundez et al.l (l2008h and the 3D radiative transfer code MC¬ 
FOST (??). Our model assumes a flared disk in hydrostatic equi¬ 
librium with a radial surface density distribution described by 
a (2(7?) oc 77 “). The best results were obtained for a physical 
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structure: Mdisk=0.01Mo, Rin =110 AU, Rout =550 AU, a sur¬ 
face density coefficient (a=1.5), and the initial abundances MO. 
Almost all modeled lines, in particular SO and CS, agree with 
the observations within a factor of 2, and the modeled prohles 
successfully reproduce the double-peak prohles from the disk 
emission. The only exception is the HCN 3—>2 that is predicted 
with an intensity around 30 times higher and the CN 2—»1 with 
a predicted intensity of 17 times higher than the observed ones, 
respectively. 

We carried out a set of modeling tasks with different initial 
abundances to study the inhuence of the initial conditions on 
our results. We took the limiting cases in which all the S is in 
CS (model MO) and all the S is in SO (model M2). We found 
that the initial conditions do not affect our result, supporting its 
robustness. The CN and HCN synthetic spectra always remain 
more intense than the observations, thus varying the initial con¬ 
ditions do not signihcantly improve our results. Only the H 2 CO 
abundance strongly depends on the initial conditions. 

Our model cannot account for the full complexity of this 
transition disk and must be considered as a guide to interpreting 
the observations. For instance, clumpyness in the disk ring would 
allow UV photons to penetrate deeper into the disk and photodis- 
sociate some molecules. A flatter disk morphology would lead 
to different chemical abundances. Rec ent observations incen- 
timeter radio emission, carried out by iRodrfguez et al.l (l2014l) . 
have shown that there is a collimated and ionized outflow in the 
inner regions of AB Aur disk. Possible shocks associated with 
this compact ionized outflow could modify the gas chemistry 
in the interaction regions. Recent ALMA observations toward 
another transition disk HD 142527 have traced HCO^ emission 
inside the ca vity, and it has been i nterpreted as accretion flows 
onto the star (ICasassus et al.llTOl^ . Finally, dust properties and 
the gas/dust ratio in transi tion disks could be differe nt from the 
standard values (see, e.g.. lvan der Marel et akllTOTih . Future in¬ 
terferometric observations would help to determine where the 
emission of different molecules are coming from and distinguish 
among all these scenarios. 
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Fig. A.I. Molecule distribution of HCN, CN, H 2 CO, HCO^, SO, and CS as a function of radius r and z/r (where z is the height over the midplane) 
for Model E. 



Fig. A.2. Distribution of H nuclei volume density {left) and temperature {right) as a function of radius r and z/r (where z is the height over the 
midplane) in Model E. 
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Fig. A.3. (Top): Model A, (Bottom): Model B. Comparison between modeled spectra and those detected by 30m telescope toward AB Aur disk. In 
colors we see the synthetic line profile obtain with our model. The blue and magenta lin es correspond to th e same model with inclination angles 
for the disk of 27° and 40° respectively. Vertical dash line indicate the Visr = 5.9 km s“* dPuvert et al.lll98^ . 
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Fig. A.4. (Top): Model C, (bottom): Model D. 
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